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PACS. 76.75.+i - Muon spin rotation and relaxation. 

PACS. 74. 70. Ad - Metals, alloys and binary compounds (including A15, Laves phases, etc.). 
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Abstract. - The evolution of the superconducting properties of the carbon-doped MgB2 
superconductors, MgF^-iCa; (x— 0.02, 0.04, 0.06) have been investigated by the transverse- 
field muon spin rotation (TF-^ + SR) technique. The low-temperature depolarisation rate, cr(0) 
at 0.6 T which is proportional to the second moment of the field distribution of the vortex 
lattice decreases monotonically with increasing electron doping and decreasing T c . In addition, 
the temperature dependence of a(T) has been analysed in terms of a two-gap model. The size 
of the two superconducting gaps decreases linearly as the carbon content increases. 



Introduction. The recent discovery of superconductivity at ^40 K in MgB2 Q has 
opened new prospects for the understanding of the microscopic origin of high-T c supercon- 
ductivity. In contrast to the complex crystal structures, multicomponent nature and compli- 
cations due to magnetism and strong electron-electron correlations of the high-T c cuprates, 
MgB2 has a simple hexagonal crystal structure (AlB2-type) comprising close-packed Mg 2+ 
layers alternating with graphite-like boron layers and sp electrons involved in the supercon- 
ducting process. Several attempts jl have been made to introduce dopants in the Mg layers in 
order to explore the relationship between T c , crystal structure and doping level. The reported 
experiments invariably lead to a decrease in T c . Successful substitution at the B sites can be 
also achieved by carbon doping || , which results in a significant contraction of the a lattice 
parameter, but affects little the interlayer separation. T c also decreases with increasing doping. 
Experimental studies of the superconducting properties of MgB2 show deviations from those 
calculated with the standard BCS theory [|],[^]. Various experiments, including scanning tun- 
neling microscopy (STM) 0, point-contact spectroscopy |Q, specific heat measurements S, 
optical H and Raman spectroscopy [|l0| suggested the existence of a secondary superconduct- 
ing gap, implying that the simple one-band approach for MgB 2 superconductivity must be 
extended. On the basis of the electronic structure, the existence of multiple gaps has been 
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invoked in order to explain the magnitude of T c in MgB2 |11|. In addition, solution of the 
full Eliashberg equations at low temperature yields different gap values for the different parts 
of the Fermi surface (~ 1.8 meV for 3D sheets and ~ 6.8 meV for 2D sheets) [^2|. Although 
there is much support for the applicability of the multiband description to MgB2, there is 
still some debate in the literature, in particular since some tunneling and NMR measurements 
show only a single gap JrJ . 

/i + SR measurements have also made significant contributions to the understanding of the 
superconducting properties of MgB2 , mainly through the precise determination of the pene- 
tration depth, A from the zero-temperature extrapolated value of the fi + spin depolarisation 
rate, er(0). Early experiments p^ j proposed that the low temperature magnetic penetration 
depth of MgB 2 shows a quadratic temperature dependence and were interpreted in terms of 
unconventional superconductivity with an energy gap that has nodes at certain points in the k 
space. However, a systematic study by Niedermayer et al. jl5| reported that the temperature 
evolution of the depolarisation rate, a in polycrystalline MgB 2 could be well interpreted in 
terms of a two-gap model. 

In the present paper we employ the transverse-field (TF) variant of the /x + SR technique 
to characterise the superconducting properties of carbon-doped MgB2 systems (MgB2- a: C :E , 
x= 0.02, 0.04, 0.06). Analysis of the experimental data in terms of a two-gap model reveals 
the effect of electron doping on the superconducting gaps and the anisotropic properties of 
the compounds. 

Experimental. - Polycrystalline MgB2- a; C a ; (x= 0.02, 0.04, 0.06) samples were prepared 
by reaction of Mg, amorphous B and carbon black at 900° C for 2 h, as described elsewhere 
The ^t + SR measurements were performed with the GPS spectrometer on the 7rM3 muon 
beamline at the Paul Scherrer Institute (PSI), Switzerland. Pressed sample pellets were 
attached with low-temperature varnish on a Ag sample-holder placed on the stick of a He 
continuous flow cryostat operating down to 1.8 K. After cooling the sample in an external 
field, H ext to temperatures below T c in order to induce a homogenous flux line lattice, positive 
muons (100% spin-polarised) with their initial spin polarisation transverse to the external field 
were implanted in the solid sample. The implanted muons come to rest at an interstitial site 
and act as highly sensitive local magnetic probes. In the presence of local fields, Bi oc , the /i + 
spin undergoes Larmor precession with frequency, uj^ = 7 M -B/ c, where 7 M /27r=13.553 kHz/G 
is the muon gyromagnetic ratio. The time evolution of the fi + spin polarisation, P/j,(t) is 
measured by monitoring the positrons preferentially emitted along the fj, + spin direction at 
the instant of muon decay. For type II superconductors, P^(t) is an oscillatory function with 
decreasing amplitude and the damping of the fi + precession signal provides a measure of 
the inhomogeneity of the magnetic field, AB in the vortex state and hence of the magnetic 
penetration depth, A. 

Results and Discussion. - For polycrystalline samples in the vortex state, the TF /z + 
spin polarisation function is approximately Gaussian (P^t) ~ exp (— ^a 2 t 2 )) and the depolar- 
isation rate, <j is proportional to the second moment of the field distribution (<r ~ 
In the case of anisotropic (A c /A a 6 >> 4) type II superconductors, a is related to the in-plane 
magnetic penetration depth, A a b via the relation ]l6| : cr[/is _1 ] = 7.086 x 10 4 A~ b 2 [nm -2 ]. 
The anisotropy of the upper critical field in single crystalline MgB2 was reported as 6 (at 5 
K) |l7j or 4.3 jig}] , while for polycrystalline samples and thin films, the reported values span 
a wide range. The above equation holds in the London limit (k = (A/£) >> 1) and in the 
absence of pinning-induced distortions of the vortex lattice. The former restriction is valid for 



MgB 2 (t; a b ~ 7 nm |19|, A„b « 95 run ITa]), while the latter is addressed below. The London 



K. Papagelis ET AL.: /^+SR study of CARBON-DOPED MgB 2 SUPERCONDUCTORS 



3 



8 



7 
6 
5 

C 4 

b 

3 



2 

1 





1 10 100 

Temperature (K) 

Fig. 1 - Temperature dependence of the fj, + spin depolarisation rate, a at H ex t~ 0.6 T for Mgf^-iCz 
(x= 0.02, 0.04, 0.06). The open circles correspond to MgBi.ggCo.o2, the triangles to MgB1.g6C0.04 
and the squares to MgB1.g4C0.06- The solid lines are fits to the two-gap model. The inset shows the 
fit of the experimental data for MgB1.gsC0.02, assuming an isotropic single gap model (dashed line) 
or a T 2 -law (dotted line). 



model predicts that the second moment of the magnetic field distribution of a perfect vortex 
lattice should be independent of H ext for A > L where L is the distance between vortices. In 
order to check the influence of pinning on the depolarisation rate, a, we have measured its 
field dependence at 5 K. Our results show that a(H) increases almost linearly up to 50 mT, 
displays a peak at ~80 mT and then decreases rapidly, reaching a plateau at higher fields 
(~ 0.6 T) in excellent agreement with the reported data for MgB 2 |L5]]. We thus performed 
our measurements at an external field of 0.6 T in order to ensure formation of an ideal vortex 
lattice and avoid underestimation of A(T). 

Fig. [l] presents the extracted temperature dependence of the TF-/i + SR depolarisation rate 
at H ext = 0.6 T for the MgB^Cz (x= 0.02, 0.04, 0.06) samples. The values of a{T) are 
derived after subtraction of the normal state temperature-independent depolarisation rate, 
oWfc (<J 2 (T) = of neas (T) - al ack ). In the case of MgBi. 98 C .o2, <r(T) increases monotonically 
as the temperature decreases below ^35 K and reaches a plateau at T < 6 K, remaining 
almost constant at lower temperatures. Analogous behaviour is observed for the other two 
compositions, except that as the doping level increases, the low-temperature plateau smooths 
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out and the extrapolated value, cr(Q) shifts to lower values. 

The experimental cr(T) dependence is reproduced well by means of a two-gap model, while 
attempts to fit the experimental data with a T 2 -law or an isotropic one-gap model led to 
unsatisfactory results (inset in fig. [I]). The two-gap model is based on the existence of two 
discrete superconducting gaps, Ai and A 2 , at T= K, both closing at T c and each associated 
with a different energy band. By assuming that the coupling between the two bands (i.e. due 
to impurity or phonon scattering) is sufficiently weak (vide infra), cr(T) can be expressed 
as 

a{T) = ct(0) - (71/7) Sa(A 1} T) - (72/7) Sa(A 2 ,T) (1) 

where 

MA,T)^ ^ f(s,T)[l-f(s,T)]de (2) 
fcs-t Jo 

f(e,T) is the Fermi distribution of quasiparticles and e the energy of the normal electrons 
relative to the Fermi energy: 

/( £ , T) = (1 + cxp(v/ £ 2 + A(Tf/k B T))- 1 (3) 

Each band is characterised by a partial Sommerfcld constant, 7; (71 +72 = 7, where 7 is the 
total Sommerfeld constant). As the Sommerfeld constant is proportional to the density-of- 
states at the Fermi level, the ratios 7^/7 determine the partial Ni(eF) for the two bands. The 
temperature dependence of the band gaps is taken from BCS theory, i.e. A(t) — A(0)5(t) 
where 5(t) is the normalised BCS gap at the reduced temperature, t — T/T c [ po| . 

The fitted parameters are summarised in Table I. For comparison the corresponding values 
for MgB 2 |l5| are also included. The evolution of the depolarisation rate er(0) with T c and of 
the superconducting gap sizes, A^O) and A 2 (0) with carbon content, x are shown in Fig. ^. 
<t(0) decreases monotonically with increasing carbon content and decreasing T c (da(0)/dT c = 
—0.25 us~ 1 K~ 1 ). The ^11% decrease in cr(0) at x= 0.06 reflects an increase of ~6% in the 
in-plane magnetic penetration depth, A a fc(0). In addition, as the carbon content increases, the 
two gap sizes at K shift to smaller values at nearly the same rate (dA^/da;= —20.4 meV, i— 
1,2). 71/7 which is a measure of the partial iV(e^) for the band with the larger gap also shows 
a tendency to increase with increasing electron doping. The superfluid plasma frequency in the 
ab plane, Wp^(= c/A a (,(0)) which is related to the charge density of the superfluid condensate 
at K is also included in Table I. 

Table I - Extracted parameters for MgB2~ x C x (x= 0.02, 0.04, 0.06). cr(0) is the low-temperature 
spin depolarisation rate, \ a b(0) is the in-plane magnetic penetration depth, Ai and A2 are the 
superconducting gaps at K, 7 (71 ) is the (partial) Sommerfeld constant, uj^ is the superfluid plasma 
frequency, n s is the superconducting carrier density and m^ b /m e is the effective mass enhancement 
of the B layers. The data for MgBi are from Ref. and the T c values from Ref. H . 
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Fig. 2 - Relationship between T c and low-temperature fi + spin depolarisation rate, <r(0) (upper panel). 
Evolution of the superconducting gap sizes, Ai and A2 at T= K with carbon content, x (lower 
panel). The closed symbols correspond to the larger gap, Ai associated with the 2D a sheets and the 
open ones to the smaller gap, A2 associated with the 3D ir sheets. 



In the London model at the clean limit, the magnetic penetration depth, A is related to 
the superconducting carrier density, n s and the effective mass, m* by 



1 4nn s e 2 
X 2 m*c 2 



1 



1 + 



(4) 



where £ is the coherence length and I the mean free path. By assuming that MgB2- K Ca; is 
in the clean- limit (£= 50 A, Z > 1000 A), the term in brackets is close to unity allowing us 
to estimate the ratio, n s /m*. The calculated values of n s /(m* b /m e ) (where (m* ab /m e ) is the 
effective mass enhancement in the B layers in units of m e ) for MgB2- x G x (Table I) are larger 
than those obtained for high-T c superconductors and decrease quasilinearly with increasing 
doping. In addition, n 3 can be converted first into a 2D carrier density, n S 2D on the B planes 
by multiplying by the interlayer distance, Co and then into an effective Fermi temperature, 
Tp through the expression, fc^Tp = fi 2 n{n S 2D / m *)- Tp decreases from —3200 K in MgB2 
to -2700 K in MgBi. 94 Co.o6 resulting in T c /Tp - 0.01. Uemura [|lj has established that 
in contrast to the conventional BCS superconductors in which T c /Tp <C0.01, the cuprates, 
fullerides and other exotic superconductors have rather high T c with respect to Tp with 
T c /Tp ~ 0.01-0.1. It is noticeable that the MgB2 superconductors also share this feature 
and are located in the Uemura plot between conventional s-wave and unconventional high-T c 
materials. 
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The electronic band structure of MgB2 has been extensively investigated |22j,|22,24|. The 
valence band of MgB2 is made up predominantly of B 2p states, which form two distinct 
sets of bands of <j(p x ,y) and 7r(p z ) type whose k dependence differs considerably. The most 
pronounced dispersion for the B p XiV states is along the T-K direction of the Brillouin zone 
(BZ), while a flat zone is formed in the k z direction (T-A), reflecting the 2D character of the 
boron lattice. These a bands are partially unoccupied creating a hole-type conduction band 
which gives rise to two 2D light-hole and heavy-hole sheets forming coaxial cylinders along 
the T-A BZ direction. The strong coupling of these holes to the optical bond stretching modes 



drives the superconductivity in MgB 2 |25|. In addition, B it states dominate at the bottom 
of the conduction band, while the Fermi surface associated with these bands consists of a 3D 
tubular network. These bands exhibit maximum dispersion along the T-A direction of the 
BZ. Theoretical calculations estimate that 44% of N(ep) comes from the 2D a cylindrical 
sheets and the rest from the 3D 7r sheets |Q. For the superconducting state of MgB2, the gap 
is nonzero everywhere on the Fermi surface and the gap values are grouped in two distinct 
regions JT2|. The larger superconducting gap, Ai on the 2D a cylindrical sheets has an average 
value of —6.8 meV, while the smaller one, A2 associated with the 3D 7r-sheets an average value 
of -1.8 meV. 

Although the electronic structure and the superconducting behaviour of MgB2 has been 
studied theoretically in detail, there is no systematic investigation concerning the effect of 
doping. In the framework of the rigid band model, it is expected |^4|] that doping of the B 
sublattice with C should lead to a shift of the Fermi level to higher energies in the region of the 
DOS minimum (pseudogap). A rigid band model estimate gives a value of —0.16 electrons/cell 
in order to fill the er bands J2j|. Electrons are also added to the it bands that lie in the same 
energy range, leading to a total doping level —0.25 electrons/cell. Thus it is expected that 
electron doping should lead to both superconducting gaps decreasing gradually. However, 
the differences in the slope — d\nT c /dV between high-pressure experiments (—0.28 A~ 3 ) on 
MgB2 and chemical substitution || (—0.37 A" 3 ), imply that the effect of carbon doping 
is more complicated due to the B layer contraction and theoretical studies beyond the rigid 
band approximation are necessary. 

Carbon substitution should also lead to an increase in the interband impurity scattering, 
which should result in the size of the a- and 7r-gaps converging to the same value. However, 
theoretical calculations |27| have shown that the particular electronic structure of MgB2 results 
in extremely weak an impurity scattering, even for low quality samples and in the presence of 
Mg vacancies, Mg-substitutional impurities and B-site substitutions by N or C. The dominant 
mechanism for impurity scattering is due to intraband scattering of the a and 7r bands with 
the scattering rate inside the 7r bands greater than that of the a bands. Intraband scattering 
does not change T c and the gap values but influences the penetration depth (Anderson's 
theorem) p8| . Our experimental data can be described well with a two-gap model in which 
impurity scattering is essentially ignored. This indicates that interband impurity scattering is 
relatively weak and both superconducting gaps are preserved, at least up to x= 0.06. Carbon 
doping also reduces both superconducting gaps with this effect relatively more pronounced for 
the smaller gap, A2 associated with the 3D sheets. The error in 71/72 (— Ai(e^) / '^(ef)) is 
quite large and hence it is difficult to extract reliably its dependence on doping. Nonetheless, 
there is a tendency for 71/72 to increase with increasing doping level. In any case, systematic 
theoretical studies are necessary to shed more light on this issue. 

Conclusions. - In conclusion, analysis of the temperature dependence of the TF-/i + SR 
depolarisation rate for MgB2- x C 2: (x— 0.02, 0.04, 0.06) shows that increasing electron doping 
is accompanied by a decrease of the low-temperature <r(0) and the spectral weight, n s /m*, 
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consistent with the decreasing density-of-states at the Fermi level. Within a two-gap model, 
both superconducting gap sizes decrease almost linearly with the smaller gap affected more 
on doping, while the interband scattering remains relatively weak up to at least x= 0.06. 
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